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Reliable determination of the complex refractive index (RI) of graphene inherently requires two independent 
measurement realizations for two independent unknowns of the real (n G ) and imaginary (k G ) components, 
i.e., RI = n G + i k G . Thus, any single set of measurement realization provides only one constraint that is 
insufficient to uniquely determine the complex RI of graphene. Tandem uses of two independent 
measurement techniques, namely the surface plasmon resonance (SPR) angle detection and the attenuated 
total reflection (ATR) intensity measurement, allow for the unique determination of the complex RI of 
CVD -synthesized graphene. The presently measured graphene RI is determined to be 2.65 + 1.27i for the 
E-field oscillating parallel to graphene at 634 nm wavelength, with variations for different numbers of L (1, 3 
and 5) remaining within ±3%. Thus, our demonstration results for the specified wavelength serve as an 
impetus to suggest the need for two independent measurement techniques in determining both the real and 
imaginary RI values for graphene. Additional efforts have been made to characterize graphene layers using 
the density function theory (DFT): this calculation provides RI G = 2.71 + 1.41i. 

The relatively significant light absorption of 2 ~ 3% per single graphene layer 1 of mere 0.335 nm implies that 
the refractive index (RI) of graphene must be complex with a significant imaginary component. Notable 
efforts (Table 1) have been made to determine the complex graphene RI using diverse experimental and 
theoretical approaches 2 " 12 , yet the available data show substantial scattering for both their real and imaginary 
components. Some of this scattering may be attributed to the graphene sample variations and/or different sample 
preparation processes, as well as to the different measurement uncertainties associated with the different mea- 
surement techniques. However, the more fundamental origins of the large data scattering are considered to be 
two-fold: (1) the graphene RI is complex- valued, having two unknowns of both real and imaginary parts at a given 
wavelength, and (2) the ultra-thin graphene layers can be naturally better characterized by using a probe focusing 
on the near-field. To uniquely determine both the real and imaginary parts of graphene RI, therefore, two 
independent sets of measurement observables will be necessary, using two discrete measurement techniques 
for the same graphene sample. Furthermore, probing near-field characterization will enhance the measurement 
sensitivity, since the surface- enhanced electromagnetic field exists enclosing the thin graphene layer. 

Employing two independent techniques requires two separate experimental systems as well as two different 
sample layouts and preparations. In contrast, notable efforts have been made by fitting experimental data from a 
single measurement technique to a dispersion model for RI of graphene. The earlier use of reflection spectro- 
scopy 2 imposed an overly simplified assumption of both real and imaginary parts of the complex RI being 
constant over the visible range, whereas its later use 3 modeled the imaginary part of the complex RI using a 
constant optical conductivity based on the earlier report 1 , yet retained an assumption of the real part of RI as a 
constant. 

Ellipsometry measures two variables (the phase and the intensity of reflected waves); however, direct inversion 
of these two measured variables into real and imaginary RI components usually requires a complex fitting process 
since all of the four variables are implicitly interconnected in Fresnel's multilayer reflectance equation. 
Spectroscopic ellipsometry 4 " 6 uses a more elaborate dispersion modeling for the spectroscopic data inversion. 
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Table 1 | Measured and calculated complex refractive index (Rl) values of various graphene samples 



Refe 



Method 



Descriptions 



Note 



Sample 



Present work 



Nietal. (2007) 2 



SPR and ATR 
measurements 



DFT 



Reflection 
spectroscopy 



Bruna & Borini (2009) 3 Reflection 

spectroscopy 

Kravets et al. (2005) 4 , Spectroscopic 
Matkovic et al. (20 1 2) 5 , ellipsometry 
Wurstbaueretal. (201 0) 6 



Jellison et al. (2007) 15 Spectroscopic 
ellipsometry 



Measurements of 9$pr ar| d 
R p /R s and simultaneous 
matching with calculation 
results 



First principles calculation by 
full potential linear muffin-tin 
orbital (FP-LMTO) within local 
density approximation (LDA) 
and generalized gradient 
approximation (GGA) 

Reflected light from a graphene 

sample by a spectrometer 
(illumination: white). 

Same as above, but the 

experimental data are from 
Ref. [1]. 

Polarization state change 
of light beam of tunable 
wavelength, when reflected 
at the graphene with 
oblique incidence 

Same as above. 



Nelson etal. (201 0) 7 , 
Weber etal. (201 0) 8 
Ye etal. (201 3) 10 

Xue etal. (201 3) 11 

Wang etal. (2008) 9 



Spectroscopic 
ellipsometry 

Spectroscopic 
ellipsometry 

Polarization- 
dependent ATR 



SPR 



Picometrology 



Klintenberg et al. (2009) 12 Density functional 
theory (DFT) 

Klintenergetal. (2009) 12 t DFT 



(2.58 ±0.38, 
1.3 ± 0.25) 

(2.63 ± 0.20, 
1.28 ± 0.10) 

(2.63 ± 0.20, 
1.28 ± 0.10) 

(2.71, 1.41) 



(2.0, 1.2) 
(3.0, 1.15) 



(2.8, 1.45) 
(2.7, 1 .35) 
(3.0, 1.35) 
(2.3,0.8) 

(2.56,2.03) 



Point-to-point fitting with no 

need of dispersion modeling 
B-spline method with no need of 

dispersion modeling 
Attenuated total internal (2.6, 1 .6) 

reflection of monochromatic, 

p- and s-polarized light at air/ 

graphene/prism with 

changing incident angle 
Reflectance of monochromatic (2.75, 0.41 ) 

and p-polarized light from 

water/graphene/ Au/ Cr/ 

glass with changing incident 

angle 

Amplitude and phase change (2.95, 1.32) 

of reflected light when a (2.98, 1 .44) 

focused monochromatic light 

traverses a graphene 

boundary. 
Same as present work (DFT) 



Same as above. 



(2.96, 1.49) 



(2.88, 1.50) 



Oscillatory behavior 

similar to that in Ref. [1 2] 



n G 



Assumed that pq and /cq 
are A-independent. 

Modeling: k G = 5.446 fim~ 

and fitted only for pq. 
Needs a dispersion modeling 
with 4 or more fitting 
parameters. 



Measurements with 
several wavelengths 



Solution optimization 

is required. 
Fitting elaboration is 

required. 
Determination of unique 

no and /cq will be 

challenging. 



Used graphene oxide 
samples with 6 fitting 
parameters to be 
determined. 

Beam profile other than 
Gaussian and slight 
defocusing should 
be noted. 

Oscillatory behavior of 
the dielectric functions 
for the visible wavelengths 



CVD, 1 L 

CVD, 3L 

CVD, 5L 

Graphene (2 
carbon atoms 
in a unit cell) 



Exfoliation, 1 L 



Exfoliation, 
Exfoliation, 1 L 



HOPG (highly 
oriented 
pyro lytic 
graphite), ZYA 
grade after 
tape cleavage 

CVD, 1 L 

Exfoliation, 1 L 
CVD, 1 1-21 



Reduced 

graphene 

oxide laid on a 

Au/Cr/glass 

substrate 
Exfoliation, 1 L 



Graphene (2 
carbon atoms 
in a unit cell) 

Graphite 



However, the reliability of the optimized result can depend on the 
stiffness of each variable and the orthogonalities between variables in 
fitting. More recent work on ellipsomery 7 ' 8 successfully demon- 
strated that assumptions on the dispersion relation modeling can 
be released, either by very precise measurements to obtain smooth 
spectral curves for the "point-by-point fitting" result, or by novel B- 
spline fitting, which doesn't require any physical assumption in spec- 
troscopic data inversion process. 

As a non-spectroscopic method, which does not require a 
dispersion modeling, picometrology 9 uses both the amplitude and 
phase change as two observables when light traverses the graphene 
edge area, but the scheme bears several error sources creating uncer- 
tainties that are not quantitatively resolved. As another non- 



spectroscopic method, the polarization dependence of optical 
absorptions under the internal reflection condition was measured 
to determine the complex RI by fitting 10 , however, the uniqueness 
of the fitted real and imaginary parts of RI is somewhat questionable 
because of the single constraint of the polarization dependence. For 
the case of a graphene flake oxide sample, fitting of the surface 
plasmon reflectance curves 11 was attempted to determine the com- 
plex RI, but again further validation will be needed for the uniqueness 
of the fitted results because of too many fitting parameters (6) for the 
single observable. 

In this report, we demonstrate successful implementation of a 
reliable and robust way to determine the complex RI of graphene 
using two independent and non-spectroscopic measurement tech- 
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Figure 1 | How to reliably measure the refractive index (RI) of graphene. (a) The experimental layout to determine the complex RI of graphene layers 
using two independent measurement realizations for the SPR angle ($spr) an d the maximum reflectance ratio R p /R s , (b) a graphene sample laid on a BK7 
glass substrate that is partially coated with 48 -nm thick Au film, (c) schematic of the two tandem measurements for 0 S pr and R p /R s under a deionized 
water environment (in lieu of air) for enhanced measurement sensitivities, (d) 0 SPR = 72.97° is measured from the maximum resonant absorption of 
surface plasmon polaritons (SPPs), or equivalently as the minimum reflectance of thep-polarized incident ray (634 nm), (e) the maximum R p /R s = 1.124 
is measured at the critical angle for total internal reflection, (f) all possible fitting solutions of (n G , k G ) that satisfy the first constraint of 0 SPR = 72.97° ', 
(g) all possible fitting solutions of (n Gy k G ) that satisfy the second constraint of maximum R p /R s = 1.124, and (h) determination of the unique 
complex RI of graphene, RI G = (w G , k G ) = 2.63 + 1.28/ from the intersection of (f) and (g), which simultaneously satisfies the two constraints. 
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Figure 2 | Measured and predicted RI G in comparison with other 
published results. The three measured RI G data respectively for 1L, 3L, and 
5L graphene samples show deviations of a mere ±3% from the averaged 
RI G = 2.65 + 1.27i. The extent of the error bars decreases with increasing L, 
primarily because of the enhanced measurement sensitivities of 0 SPR with 
thicker graphene samples while the measurement uncertainty for Rp/R s 
remains unchanged. The predicted result of RI G = 2.71 + 1.41/ (*) using 
the density function theory (DFT) agrees fairly well with all of the measured 
RI G data. RI G data from previous efforts include the reflection 
spectroscopy 2 ' 3 , the spectroscopic ellipsometry for graphene 4-6 and for highly 
oriented pyrolytic graphite 151 , more recent ellipsometry that alleviate the 
need of dispersion modeling 7 ' 8 , the picometrology 9 , absorption under 
internal reflection 10 , surface plasmon reflectance for reduced graphene oxide 
(rGO) 11 *, DFT calculations for graphene 12 and for graphite 121 . More detailed 
descriptions on these published data are shown in Table 1, and further 
elaboration on their uncertainties is presented in the main text. ( ': Results 
for graphite samples; *: Results for reduced graphene oxide sample; All other 
data: Results for CVD -synthesized or exfoliated graphene samples). 



niques: (1) detection of the maximum absorption of the p-polarized 
incidence at surface plasmon resonance (SPR) angle 0 SPR by surface 
plasmon polaritons (SPPs) of a thin Au layer coated with graphene, 
and (2) detection of the reflectance ratio of R p /R s at the critical angle 
when the incident light is attenuated by the graphene sample without 
an Au interlayer. These two independent observables allow for the 
unique determination of the real and imaginary parts of the complex 
RI of graphene. The main point that we would like to report here is 
that the proposed tandem technique, unlike any previously pub- 
lished techniques, can provide two independently measured vari- 
ables, allowing the straightforward determination of both the real 
and imaginary components of graphene RI with no fitting 
elaboration. 

Results 

The main design of the experimental system is comprised of two 
rotatable arms and a prism assembly placed near the pivot point 
(Fig. la): (1) the left arm consists of the incident LED light source 
of 634 nm wavelength with a FWHM (full width half maximum) of 
22 nm, a polarizer and collimating optics, (2) the right arm consists 
of an EMCCD (electron-multiplied charge-coupled device) camera 
and microscopic lens assembly, and (3) the top surface of the BK7 
prism accommodates an index matching BK7 glass slide that is par- 
tially coated with a 48-nm thick gold layer (Fig. lb). The gold coated 
surface is used to measure 0 S pr and the uncoated BK7 glass surface is 
used to determine the maximum R p /R s (Fig. lc). Note that the water 
environment is provided for the graphene sample to enhance the 
angular measurement sensitivity by increasing the magnitude of 
the SPR angle ranges, in comparison to an air environment 13 . 

When the surface plasmon polaritons (SPPs) - the surface charge 
density waves of the p-polarized incident light - are resonantly 
coupled with the free electrons in the Au layer, the resulting reflec- 
tance R p falls to a minimum 14 at 0 S pr (Fig- Id)- For the four-layered 
structure (BK7 glass-48 nm Au layer-graphene layer- water), 9 SPR is 
given as a function of both the real (n G ) and imaginary (k G ) compo- 
nents of RI of graphene (Fig. If) based on the modified Fresnel 
equations (see the first Supplementary note online). For the case of 
3L, for instance, the experimentally measured 0 S pr = 72.97° satisfies 
all combinations of (n G , k G ) as marked in blue in Fig. If, and thus, it is 
clear that the determination of a unique set of n G and k G requires an 
additional constraint from another independently measured set of 
data. 

The additional constraint comes from the measurements of the 
maximum R p /R s as schematically illustrated in Fig. lc. When an 
absorbing layer such as graphene disjoins the water-glass (BK7) 
interface, which is subjected to total internal reflection, the evan- 
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Figure 3 | Uncertainties associated with the single constraint of R p /R s data, (a) The measured maximum Rp/R s , i.e., Rp/R s * = 1. 124, is shown as the solid 
curve (black) on the calculated 3-D contour (see second Supplementary note online), and its projection onto the (n G , k G ) plane (the gray curve) 
corresponds to all of the fitted solution candidates for RI G = n G + ik Gy under the given single constraint, and (b) three arbitrarily selected pairs of (n G > k G ) 
among these solution candidates predict an identical R p IR s curve (see second Supplementary note online) and all of them fit to the experimental R p IR s 
data, which demonstrates the difficulty in uniquely determining the complex RI of graphene using the R p IR s data alone. 
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Figure 4 | Uncertainties associated with the single constraint of # S pr data, (a) The measured SPR angle, i.e., 0 S pr = 72.97°, is shown as the solid curve 
(black) on the calculated 3-D contour (SI1) and its projection onto the (n Gy k G ) plane (the gray curve) corresponds to all of the fitted solution candidates 
for RI G = n G + ik G , under the given single constraint, and (b) three arbitrarily selected pairs of (n G , k G ) among these solution candidates predict an 
identical surface plasmon reflectance curve (SI1), thereby an identical $spr> an d all of them fit to the experimental 0 SPR data, which demonstrates the 
difficulty in uniquely determining the complex RI of graphene using 0 SPR data alone. 



escent wave is partially attenuated and the resulting reflectance 
reduces from unity. Furthermore, the amount of reduction in reflec- 
tance is different for thep-polarized and the s-polarized incident light 
and the ratio of R p /R s is dependent upon the incident ray angle 
(Fig. le). The maximum R p /R s at the critical angle is selected as a 
second independent observable to provide an additional constraint 
for the unique determination of (rc G , k G ). As expected, the measured 
value of maximum R p /R s = 1.124 satisfies multiple pairs of (n G , k G ) 
as shown in yellow in Fig. lg, and again, unique determination of 
complex RI is not possible by this constraint alone. Now, overlapping 
Fig. If and Fig. lg creates the intersection region that simultaneously 
satisfies both constraints of 0 SPR and maximum R p /R s (Fig. lh), and 
therefore, the complex graphene RI is determined to be RI G — 2.63 + 
1.28/, for the case of the three-layered graphene (L = 3), which is 
indeed the average of the two matched pairs corresponding to the 
intersection. 

The resulting RI G data are shown in Fig. 2 as n G + i k G = 2.58 + 
1.30/, 2.63 + 1.28/, and 2.73 + 1.23 i for 1=1,3, and 5, respectively. 
The average of these three gives RI G = 2.65 + 1.27 i and the varia- 
tions for different L remain within ±3% from the average. The 
smaller uncertainty range for the higher L is attributed mainly to 
the higher RI sensitivity of 0 SPR with increasing L, and this will be 
further elaborated in the Discussion section. Calculations are also 
conducted to predict RI G using the density function theory (DFT) 
and the present result of RI G = 2.71 + 1.41i is marked with the star 
symbol in Fig. 2. In comparison with the published DFT results 
of 2.96 + 1.49/ for graphene 12 and 2.88 + 1.50/ for graphite 12t , 



the discrepancies are attributed to the minor differences of the 
calculation details such as the cutoff energy and the atomic potential 
modeling, while the dielectric functions agree well for the photon 
energy up to 20 eV. 

Also presented in Fig. 2 are published graphene RI data obtained 
from one of the single measurement techniques imposed with addi- 
tional conditions and/or assumptions. These data are scattered for 
the real part, 2.0 < n G < 3.2 (±23% scattering), as well as for the 
imaginary part, 0.78 < k G < 1.6 (±34% scattering). The two results 
from reflection spectroscopy 2,3 show large deviations in their real 
parts of RI G in particular. The former used a constant graphene RI 
model 2 that was later proven to be inappropriate 4,5 ' 7 " 912 . The latter 3 
used an incomplete dispersion model under the assumption of the 
universal optical conductance 1 and assumed a constant of the real 
part of the graphene RI as n G = 3 for the best fit. Later, this was found 
to be inaccurate from the more up-to-date DFT analysis as well as 
from some experimental findings 4,5,7 " 9,12 . The deviation of the ellip- 
sometry result 6 from those using a similar ellipsometry technique 4,5,7,8 
may be attributed to the simplified Drude model used for the data 
inversion, whereas the latter papers used the Cauchy model 4 , Fano 
resonance model 5 , point-to-point fitting/Lorentz oscillator model 7 , 
and B- spline method 8 , respectively. 

Also, the discrepancies of some experimental results may be due to 
their uses of different samples other than pristine graphene layers, 
such as reduced graphene oxide (rGO) flakes 11 * or highly oriented 
pyrolytic graphite (HOPG) 15f . Graphene oxide flakes are subjected to 
form segmented layers with impurities that can lead to reduction of 



Rs 5L 





Figure 5 | Measurement sensitivity of # SPR vs. the uncertainties of (n G , k G ). The predicted 0 SPR contours ( SI 1) with consistent 0.2° increments are shown 
for L = 1,3, and 5, and the slope of the contours above the (n G , k G ) plane increases with increasing L. This implies that, given a fixed measurement 
sensitivity of ±0.1° for 0 SPR , the corresponding uncertainties in determining (n Gy k G ) are reduced with increasing L. 
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3.55 

Figure 6 | Orthogonality of the two independent observables, # SPR and R p /R s *. (a) The orthogonality of the two sets of solution pairs (n G , k G ) from the 
two constraints, 6 SPR and £ p /£ s *, is defined as the sine function of their intersection angle a on the projected plane, and (b) the maximum orthogonality of 
unity for a = 90° provides the strong independence of the two constraints ensuring the unique determination of the complex RI of graphene. All of the 
three measured RI values for L = 1,3, and 5 fall on the desirable maximum orthogonality of unity. 



electron mobility, and equivalently, reduced imaginary part of RI 
(k G ). In contrast, the agility of electrons in the highest-grade 
HOPG can result in high optical conductivity and a noticeably large 
k G . 

Discussion 

Figure 3 a illustrates the fitting uncertainties of the complex graphene 
RI to the measured R p /R s for the case of the 1L graphene sample. The 
three curves correspond to three arbitrarily selected pairs of (n G , k G ) 
from the solution candidate pool (Fig. lg) and all three curves fit to 
the measured maximum R p IR s * = 1.124. This implies that any of the 
three pairs, (3.45 + 1.00/), (2.63 + 1.28/) or (1.90 + 1.90/) can be 
accepted as a fitted value for the complex graphene RI. Indeed, pro- 
jection of the measured constraint R p /R s * = 1.124 onto the n G -k G 
plane (Fig. 3b) shows infinite pairs of (n G , k G ) that satisfy the mea- 
sured constraint. Measurements of SPR reflectance (Fig. 4a) also 
show that all three curves generated with three arbitrarily selected 
(2.80 + 1.85/), (2.63 + 1.28/) or (2.58 + 0.70/) fit closely to the 
measured and 0 SPR = 72.97°. Again, the n G -k G projection of the 
corresponding contour of measured 0 SPR = 72.97° provides only a 
partial constraint that any of these pairs of (n G , k G ) located on the 
projection band can be accepted for the fitting solution. 

The contacts between graphene and a variety of metals were dis- 
cussed by theoretically investigating the Fermi level shift in the con- 
tacted graphene from the freestanding one, and calculating the 
electronic structure and electrostatic potential 16 . Their findings 
include that for the Au-graphene contact, there lies a potential bar- 
rier between them, which will cause extra contact resistance, result- 
ing in only weak physical contact. Thus, we believe that graphene 
maintains weak physical contact with the gold substrate and that the 
opto- electric properties of graphene, including its complex RI, are 
not altered by the contacting Au substrate. Furthermore, we have 
found that a number of research groups have used different sub- 
strates contacting graphene, assuming no significant changes in gra- 
phene's RI properties for the case of amorphous quartz substrate 4 , 
GaAs substrate 6 , poly-dimethyl-siloxane (PDMS) 10 , and also for Au 
substrate 11 . 

For AA-stacking of multilayer graphene, i.e., atom-to-atom 
arrangement of graphene layers, no measurable distinctions in RI 
were observed with the number of graphene layers 17 . For twisted 
bilayer graphene including AB-stacking 1819 , in contrast, the 
reflection contrast of multilayers was found to slightly vary with 
the existence of the 2 nd graphene layer; however, more consistent 



and quantitative conclusions are still under examination. Note that 
these findings were relevant for mechanically exfoliated graphene 
layers from highly ordered pyrolytic graphite (HOPG), whereas 
the present study examined CVD- synthesized graphene layers that 
were created by the crystallized growth of many randomly oriented 

a 30 




Energe (eV) 

Figure 7 | Ab-initio DFT calculation results for the complex dielectric 
function of graphene. The real (a) and imaginary (b) parts of the complex 
dielectric function of pristine graphene are calculated using the density 
functional theory (DFT). The dielectric function s xx is for the electric field 
oscillating parallel to the graphene layer and e zz is for the electric field 
perpendicular to the graphene layer. Note that our proposed scheme, like 
all other available experimental methods, is able to conduct measurements 
only for the xx-components of the complex RI where the electric field 
oscillates in parallel with the graphene layer and is unable to detect the zz- 
components. For the incident wave of 634 nm corresponding to 1.96 eV, 
the real and imaginary parts of the relevant dielectric function are given by 

= 5.36 + 7.64/ and the resulting complex RI of graphene is given by 
^= n G + ik G = 2.71 + 1.41/. 
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submicron-sized grains. Therefore, each layer of CVD- synthesized 
graphene is considered to be randomly oriented, and the multi- 
layered graphene obtained by repeated transfer of each layer does 
not provide any consistent stacking orientations. 

The increased total thickness of graphene with increasing L creates 
steeper contour surfaces above the n G -k G plane (Fig. 5), which in turn 
contributes to lowering the measurement uncertainties. The smaller 
error bars of the measured complex RI data with increasing L, as 
shown in Fig. 2, are attributed to this enhanced sensitivity with 
increased total thickness of the graphene samples. It is worthwhile 
noting that the multi-layered graphene samples are prepared by 
physical stacking of multiple CVD -synthesized layers, one at a time, 
through multi-step processing such as spin coating, baking, and 
etching of PMMA and use of chemicals 20 and extra care must be 
taken to minimize left-over impurities and contaminants. It is worth 
noting that a small amount of impurities would make discernible 
changes in the measured 0 S pr resulting in an overestimated n G since 
the SPR is known as the most sensitive detection tool for the real 
component of RI 21 . In addition, the Raman signal can provide quite 
reliable information on imperfections, disorders, and grain bound- 
aries on graphene. However, the PMMA residues that may result 
from incomplete washing after the transfer process are not readily 
detectable by Raman signal since they are Raman inactive 22 . 

Independence of the two observables, 0 SPR for minimum SPR 
reflectance and the maximum R p /R s for attenuated reflection, can 
be assessed by their mutual orthogonality on the n G -k G plane 
(Fig. 6), which is defined as a sine function of the intersection angle 
of the two projected contour bands. Thus, the orthogonalities of near 
unity of the present cases of L = 1, 3, or 5 supports fairly unique 
determinations of the graphene RI within acceptable experimental 
uncertainties. Also, the fitting uncertainties for both SPR angles and 
total reflectance ratios associated with FWHM of 22 nm are shown 
to be less than 0.3% deviations from those with zero FWHM, based 
on the reflection calculations presented in Supplementary 
Information. 

Methods 

Preparation of graphene sample on a partially Au-coated glass substrate. Single 
graphene layers were grown on thin Cu foils (25 |im thick, 99.9999% purity by Alfa 
Aesar Inc.) using a commercially available low-pressure-high-temperature (LPHT) 
chemical vapor deposition (CVD) system by ScienTech Inc. A portion of the slide 
glass substrate was masked with a cover slip using Kapton tape (DuPont and Kapton 
Inc.) before the adhesion layer of Cr (4 nm) and Au thin layer (48 nm) were 
sequentially sputtered. Then, the cover slip was removed to recover the region of 
uncoated glass substrate. Finally, CVD graphene layers were transferred onto the 
substrate using a conventional PMMA (poly-methyl methacrylate) process 20 and the 
final test sample was completed as shown in Fig. lb. 

Incident angle measurements. The angle of incidence was measured by a Pro 3600, 
Mitutoyo digital protractor with 0.01° minimum reading scale. The tilt angles of the 
two arms of the experimental system (Fig. la) were measured with respect to the 
horizon and the actual incident angle to the graphene sample placed on top of the BK7 
prism was accounted for the refraction of the incident ray into the prism. The 
measurement increments ranged from 0.1 to 0.7° to accommodate minimum 
distinguishable measurements of the reflectance data. 

Data acquisition and recording. The reflected light intensity was detected by a 
Princeton Instrument ProEM 512b EM-CCD (electron-multiplied charge-coupled 
device) with 512 X 512 16-micron-square pixels at a 16-bit data transfer rate. The 
non-linearity of the CCD due to the low- noise was kept under 1% and the coolant 
temperature was kept at — 70°C. 

Calculations of complex RI of graphene using Density Function Theory (DFT). 

Our calculations for complex graphene RI used the Vienna Ab-initio Simulation 
Package (VASP) to determine the dielectric constants for both the in-plane (xx) and 
cross-plane (zz) directions (Fig. 7). Then, the complex RI value of graphene can be 
readily obtained by taking a square root of a complex dielectric constant at a specific 
energy (equivalently, wavelength" 1 ). The charge density distributions were first 
calculated for use in the electronic band structure and dielectric function calculations. 
In calculating the charge distributions, graphene was set to have a constant inter- 
atomic distance of 1.42 A and the tetrahedron method with Blochl corrections was 
used for partial occupancy with broadening of 0. 1 eV. The projector-augmented wave 
(PAW) method was adopted with ultra-soft pseudo-potentials and energy cutoff of 



400 eV. The K-point grid was generated by automatic generation of T-centered 36 X 
36 X 12 mesh and the electronic band structures were calculated with Gaussian 
smearing along the high symmetry points of the irreducible Brillouin zone of 
graphene, T-M-K-r, with 20 segments in each line. The present result shown by the 
star symbol in Fig. 2 comes from consideration of the dielectric constant at 1.96 eV, 
which corresponds to the center wavelength of 634 nm of the LED light that was used 
for the experiment. 
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